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Abstract 

For systems of one-component interacting oscillators on the d-dimensional lattice, d > 
1, whose potential energy besides a large nearest-neighbour (n-n) ferromagnetic translation- 
invariant quadratic term contains small non-nearest-neighbour translation invariant term, an 
existence of a ferromagnetic long-range order for two valued lattice spins, equal to a sign of 
oscillator variables, is established for sufficiently large magnitude g of the n-n interaction with 
the help of the Peierls type contour bound. The Ruelle superstability bound is used for a 
derivation of the contour bound. 



1 Introduction and main result 

Let's consider the system of one- dimensional oscillators on the d- dimensional lattice Z"^, with 
the potential energy (on a set A with the finite cardinality |A|) 

U{qA) = T.Hq.)-2dgql)+g (Q^ - Qy? + U' (qA) , g>l (1.1) 

xgA \x—y\=l,x,y^A 

Here q^ the oscillator coordinate taking value in M, qx =(qx,x G X), the one-particle 
potential (external field) m is a bounded below even polynomial having a degree degu = 2n, U' 
is an even translation invariant function such that U satisfies the superstability and regularity 
conditions, |x| is the Euclidean norm of the integer valued vector x, c? > 1. 

Let <>A, <> denote the Gibbs classical or quantum average for the system confined to A 
and the system in the thermodynamic limit, i.e. A = Z'*, respectively. 

For classical systems 

< Fx >A= / Fxiqx)e-''''^''^UqA, Za = J e-^''^'^^^ dq^. 

If Fx is the operator of a multiplication by the function Fx{qx) then the quantum average 
is given by 

< Fx >A= Z^'TriFxe-^"''), Zj, = Tr(e-'^^') 

where = 9^ + U{q\), and is the partial derivative in q^. 

The corner stone of proving an existence of Iro, using generalized Peierls argument, is the 
following contour bound 

< n X^X.- >A<e-^''', (1-2) 

<a;,a;'>er 

where F is a set of nearest neighbours, |F| is the number of them in it. 
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Xt = X{0,oo){qx), Xx = X(-oo,0)(ga:), 

X{a,b) is the characteristic function of the open interval (a, b). 

The bound (1.2) was earher derived in [BF],[FL] for several classes of classical ferromagnetic 
systems or classical systems with the nearest-neighbour pair interaction (see,also,[Si],[SH],[BW]) 
If one puts Sx = signqx, then taking into account that xt^~^ — ![-'- + one obtains 

4 < XxXy >A= 1+ < Sx >A - < Sy >A - < SxSy >A . 

Since the systems are invariant under the transformation of changing signs of the oscillator 
variables we have 

< SxSy >A= 1 - 4 < xtXy >A • 

Now in order to prove the ferromagnetic long-range order for the spins Sx one has to show that 
the average in the rhs in the equality is strictly less than j. This can be proved with the aid of 
the following lemma ([GJS], [FL]) 

LEMMA 1.1 

If the bound (1.2) holds, d > 1, and is sufficiently small then there exist positive 
numbers a, a' such that 

< xtXy >< a'e-"^. (1.3) 

So, if one shows that E can be made arbitrary large while increasing g or /3, then the Iro 
for the above spins will be proved. 

In this paper we prove the ferromagnetic Iro for the systems, in which interaction is neither 
ferromagnetic nor n-n, but essentially ferromagnetic for sufficiently large g (see Remark 5). We 
establish (1.2) for the simplest polynomial u{q) = riq"^"- with the help of the Ruelle superstability 
bound [RI] and show that E in (1,2) is positive and growing for growing g, or more precisely 

£; = eo - 2-Mn(87r-%) - £;o, eo = [^"2d(r/n)-^]^, (1.4) 

where Eq depends on g, (3 and is found from the superstability bound (for the rescaled and 
translated correlation functions). Eq is bounded in g for classical systems and grows slowly for 
quantum systems (see Lemma 1.2). 

The proposed technique is based on the precise knowledge how the constant, defining Eq 
in the superstability bound, depends on the potential energy (Theorem 2.1). It is inspired by 
the technique proposed in [AKR] for quantum ferromagnetic systems, which by rescaling of the 
oscillator variables, can be reduced to the above systems with the pair quadratic infinite-range 
interaction 

U' ^ Cx-y{qx - qyf, u{q) = r]q^. 

x,yeA 

In this paper a complicated version of (1.2) is proposed and the Iro is proved for Gibbs loop 
path system associated to the quantum system via FK formula and unit spin which are signs of 
the averaged Wiener path. A small parameter, appearing in the potential energy, determining 
a depth of the symmetric wells of the external potential is not associated with the magnitude 
of n-n interaction in it. 
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Our approach stresses the necessity to consider a large magnitude of n-n interaction which 
determines the depth of the symmetric wells of the external potential u{q) — 2dgq^ (cq, — cq are 
the only its real minima). 

Proofs of an existence of an order parameter for ferromagnetic quantum oscillator systems 
with n-n interaction, which are based on the reflection positivity, can be found in [KP],[BK] 
(see also [DLS]). Vanishing of the order parameter in the quantum limit (mass is vanishing) is 
estabhshed in [VZ]. 

THEOREM 1.1 

Let the potential energy of the one-component oscillator classical or quantum system is 
given by (1.1), u{q) = rjq'^"'. Let, also, U' be a translation invariant and an even function such 
that the condition of superstability and regularity hold for it 

U'{qA)>-J2[B/\qx\' + B'], 2</<2n, 
\W'{qxr,<i.J\ = \U'{qx,ux,) - U\qx,) - U'{qx,)\ < 

x£Xi,yeX2 

where /' < |, / > 2; /' < |, / = 2 for non-negative U' in the second inequality and /' < in both 
inequalities if U' is non-positive; B,B',^\x\ are non-negative constants, ||^'||i = S ^|a;| < °o 

and the summation is performed over Z*^. 

Then there is the feromagnetic Iro in classical and quantum systems for the spins for 
sufficiently large g : g » 1, i.e < SxSy » 0. 

Since Sx are scale invariant and their average is not changed after rescaling of oscillator 
variables, we can deal with the rescaled by g~2 variables and the potential energy Ug 

UgiQA) ^ J2 + E {qx-qyf + U'{g-kA), (1-5) 

xeA \x-y\=l,x,yeA 

where 

u^iq) = - 2dq'. 

The correlation functions or reduced density matrices generated by Ug will be denoted by 

Pg- 

The main idea of the proof originates from the inequality 

< n X^Xr'>A<(e'eo)^e-°l^l<e«-->A, (1.6) 
<x,x'>er 

where e' is a positive constant, Cq is a growing function of g, the expectation value is determined 
by Pg and 

1 4 

QgAQA)= E QgilxAy), Q9{<lx,qy) = —{{.qx-qx'f + i;{\Qi-el\ + \ql,-el\)}. 

<x,x'>&T ^0 

Here we used the inequality 

X^{(lx)x~{<lx') < {e'eo)h-^° exp{Qg{qx, qy)]}. (1.7) 
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Theorem 1.1 will be proved if we prove the following lemma. 



LEMMMA 1.2 

Let the conditions of Theorem 1.1 be satisfied. Let, also, eo be given by (1.4). Then there 
exists a function Eq^q) on the interval [1, oo) such that 

< e^"'^ >< el^l^°. (1.8) 

For the classical systems Eq is a bounded function on the interval [1, oo). 
For the quantum systems if 



k(g) = (1 + e-VlP)-\l - e-VI/3) _ ^eo-\/^ > (1.9) 

o V 71% 



20 

y 

then there exists a bounded continuous functions E^{g) on [1, oo) such that 

Lemmas 1.1, 1.2, i.e. (1.3) and (1.10) prove Theorem 1.1 since eo grows faster than ^/g. 
Function E^ in the Lemma is defined by (3. 9-11), (3. 14). 

In the second and third sections we'll give the proof this lemma for classical and quantum 
systems, respectively. Proofs of Lemma 1.1 and (1.7) are standard and will not be give here 
(see [GJS],[FL],[AKR]). 



2 Lemma 1.2 via superstability argument. Classical sys- 
tems. 

For classical systems with the rescaled potential energy 

< Fx >A= ZX' I Fx{qx)e-^'''^'^^^dqA = J Fx{qx)Pg{qx)dqx, 

p^iqx) = Zj,' J e-^^«(«-)d^A\x, Z^ = J e'^'^^^^-^dqA. 

Here the integration is performed over R\^\ and are the correlation functions. By pg we'll 
denote the correlation functions in the thermodynamic limit. 

Changing the varibles qx ^ qx ~ ^o, in the integral in the right-hand-side of (1.6) and using 
the translation invariance of the Lebesque measure we obtain 

< e'^s-r >= j pg{qY + eo) exp{(5g,r(gr + eo)}(igr, qv = {qx, qy] <x,y>eV), (2.1) 

10 8 

Qfl,r(gr + eo)< {^(^^ + + o (l^^l + gx + eo = (?x + Cq, x e X). 

The polynomial Q becomes bounded in g if it is translated by Cq. As a result, we have to 
prove that the correlation functions, translated by eo, in the limit of growing g satisfy the usual 
superstability bound. 

It is not difficult to check that if eo is given by (1.4) then 
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Prom this we immediatly see that 

^ . 

^0(^ + 60) =pg(g) + - 6', b=2dn~U2n(n-l)-n), b' = 2d 

where is a bounded below polynomial in eg ^ and q (the linear term proportional to eo is 
absent in it) 

Pg{q) = 2dn 'Y.— r^l 4 ■ 

s=3 

Now we have to establish the accurate superstability and regularity conditions for the trans- 
lated by eo potential energy. 

The superstability bound is given by 

Ug{qx + eo) > E - l^lBg, Bg^b' + B' , (2.2) 

xex 

where 

u,{q)^{ul{q + eo) + b')-Bg-'^^'\q\K 

For a non-negative U' 

Ug{q) =wJ(g-Feo) + b'. 

Let's put 

U„g{qx) = Ug{qx + eo) - X! + \^\Bg, (2.3) 

u*g = Ug-Vg, Vg{q) = q^ + g~'^^^'\q\K 

Eg diverges if g tends to infinity since b' . We can add |A|i?g to the potential energy since the 
expression for the correlation functions is not changed after this. 
Then the following superstability condition holds 

U.{qx) > E '"M- (2.4) 
xex 

The regularity condition, also, holds 

\W^g{<lXi]qX2)\ = \U*g{qXiUX2) - U^giQXi) - t^*g(fe)| < 

<l E '^\x-y\[Vg{qx)+Vg{qy)], H X2 = 0, (2.5) 

^ xeXi,yeX2 

where "^i^i = 25|^|,i + ^'j^i. 

Applying |X| — 1 times the regularity condition the following important condition is also 
derived 

U.g{qx) < E Ug{qx), Ug{q) = U.g{q) + miiVgiq)- (2.6) 
xex 

From the definition of the functions determining Ug, taking into account that Ug{q) — Ug{q), 
we derive 
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ilgiq) ^B'+{1+ \Mi)q^ + {1 + B+ 

Let's put 

ptgiQA) = exp{/3 J2 u*g{Qx)}Pg{qA + eo). 

Then p^g are expressed in terms of after adding to Ug the large in g terms independent 
of oscillator variables 

P%{qx) = Z-j^ I e-^^-(^-)/x.,(dgA\x), = / e-^^-^'^-^/.^lrfgA), (2.7) 

where 

IJ'*g{dqY) = exp{-^ u^gi(lx)}dqY- 

As a result of the superst ability and regularity conditions for U^g the following theorem is 
true [R]. 

THEOREM 2.1 

Let the condition (2.4-5) hold for a positive polynomial Vg and the function u^g be such that 
the measure is finite. Then for arbitrary 0<3e<l,r>0 for the correlation functions 
defined by (2.7) the following (superstability) bound is valid 

ptg{qx) < exp{- E [/3(1 - ^e)vg{q,) - Co(/-i^, 4,J]}, (2.8) 

where c is a positive continuous monotonous growing at infinity function, 

\q\<r 

lu^ j exp{-/?[(l - ^e)vg{q) + u{q)]}dq. 

We formulated the Ruelle result in a more explicit form in order to trace the dependence in 
g in all the terms. 

(2.1) and theorem 2.1 yields 

< eQ«.r >< giri^o^ E^^E' + eM. e^g) ^ co{l;l^JuJ, (2.9) 

E'^2\nl exp{-/?(l - 3e)vgiq) - (3u,g{q) + ^ ^_\q\'^dq. 

As a result, (1.2) holds with E given by (1.4). From the conditions of the Theorem 1.1 it 
follows that and exist in the limit of vanishing g~^. Here we have to rely on the following 
significant equalities 

lim iu^Aq + Co) + b') = bq^, lim v„iq) = kq^, b > Ad, 

where k = 1 or k = 2. From the inequalities (|g| < r) 

Ug{q) < B' + (l + ||*||i)r' + (l + B + ||*||i)y-^+''r', (2.10) 
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u.g{q) < Ug{q) + Vg{q) < Pg{r) + + Vg{r) + Bg-^+^'r^ (2.11) 

it follows that 

te,)-' < r-^e^^^'^'^^ (2.12) 

p+(r) = pg{r) + + (2 + ||*||i)r' + (2 + B + H^HO^-^+'V. (2.13) 
Polynomial />^(r) is uniformely bounded in g 

Pg{r) < p°(r) + B' + {2 + ||*||i)r^ + (2 + S + 1 1*| |i)r^ = p(r), (2.14) 

0/ ^ _i^s!(2n-s)! , 
p"{r) = 2dn ^ 2^ — j — -r . 

s=3 

So, e^{g) is a bounded function. 

Classical part of Lemma (1.2) is proved. Application of Lemma 1.1 completes the proof of 
Theorem 1.1 for classical systems. 

3 Lemma 1.2 via superstability bound. Quantum sys- 
tems. 

For quantum rescaled systems the Gibbs average of the operator Fx of multiplication by the 
the function Fx{qx) is determined by the reduced density matrices (RDMs) Pg{qx\qx), 

< Fx >A= ZX'TriFxe-''''^) = 
= Zl\^P j Fxiqx)e-^''HqA;qA)dqA = J Fx{qx)p'^{qx\qx)dqx^ (3.1) 

p''g{qx\qx) = (V^)l^l / p'',MPf,,g,{dwx\ p'^M = Z^^ j e'^"^--) P,{du^\x). (3.2) 

where uo = {q,w) ^ Q* = RxQ, Q is the probability space of Wiener paths(i(7 G fi), Pqq,{dw) 
is the Wiener (conditional)measure concentrated on paths, starting from q and arriving in q' 
at the time t, Po{duj) = y^dqP^'^ldw) , 

3/3 



Ug{ujA)=g-^ j Ug{wA{t))dt = J Ug{wA{gt))dt 



In deriving the formulas we applied the Feymann-Kac formula to the kernel 
e~^^ {Vg^qx', Vg'^q'x) of operator e~^^ and the relation 



(3.3) 

which follows from 
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eMtd'K^'q; y/FW) = (47r0-^ exp{-^^} = ^eMt9d'}{q; q'). 
The rescaled Hamiltonian is given by 
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x€A 

In order to prove Lemma 1.2 one has to estimate Pg{qx + ^ollx + Cq)- 
From the translation invariance of the conditional Wiener measure and the measure Pq it 
follows that 

P^iqx + eo\qx + eo) = (V^)"" / Pgi^x + eo)Pf,,,,{dwx), 

where cu + Cq — w{t) + Cq, w{0) — q,t & [0, t\. 
As a result 

< e«-^ >A= {^f\ J e«-^(^^+^°)pJ(a;x + eo)dqrP^^^^g^{dwr), 

It is evident that Lemma 1.2 can be proved now with the help of the analogue of the supersta- 
bility bound for pgiojp^ + cq) which was proved by Park [P]. 

In order to prove the analogue of Theorem 2,1 one has to derive the sup erst ability and 
regularity conditions for U^g{ujjs) = Ug{uA + eo). But now it is easy since in the previous section 
we estalished them for Ug{qA + eo). 

So, let by u^-.g{uj), Vg{uj), Ugiyj) be denoted the corresponding functions, depending on w{gt), 
being integrated by dt on the interval [0, (3] and 

(ua) ^ J U^g{wA{gt))dt. 



where U*g{w\{gt)) is defined by (2.3)(instead of lu w may be written). Then 

U^LVx) > E ''aM- (3-4) 
xex 

\W^g{toxi;uJx2)\ = \U*g{u;xiuX2) - U^g{uxi) - C^*g(?X2)| < 

<\ E ^k-2/iKK) + ^9K)]> XinX2 = 0, (3.5) 

^ x&Xi,yeX2 

U^{uJx) < E Ug{uJa:). 

xex 

Let's put 

ptgM = explE u^gMjPgi^^A + eo). 
xeA 

Hence 

piicox) = Z;l J e-^-(-^)p,o(c?^A\x), ^*A = / e-^-(-^)p,o(c^c^A), (3.6) 

where 

P^oiduy) = exp{- E u*g{oUx)}Po{duY). 
As a result of (3.3-5) the following theorem is true. 
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THEOREM 3.1 

Let the condition (3.4-5) hold for a positive polynomial HgiQ) and the functional u^,g be such 
that the measure P*o is finite. Then for arbitrary 0<3£< l,r>0 for the correlation functions 
defined by (3.6) the following (superstability) bound is valid 

p^iux) < exp{- Y: [(1 - ^e)vgM - coil^l^, luj)]}, (3.7) 
where is a positive continuous monotonous growing at infinity function, 

Ir,u = e-5^ll*lli^«''^7o«, Vg^r = ess sup Vg{cu), 

hu.g = J ex.p{-[Ug{uj) + u^g{u;)]}Po{du;), = j exp{-[(l - 3e)vg{u;) + u^g{u;)]} Po{du;) , 

Q* 

r 

where ft; ^ {lo E fl* : \w{t)\ < r}. 

The proof of this theorem does not essentially differs from the proof of Theorem 2.1. 
Proof of Lemma 2.1. 

Theorem 3.1 yields the following equahties 

< e^-r >< giriEo^ Eo^E' + e^g), 6.(5) = co(/-i^, 4 J, (3.8) 

E° = 21n exp{-/?[(l - 3e)vg{w) - u,g{w)] + ^^q' + ^\q\}Pf^{dw)dq. 

e^,{g) is a bounded in g since Ug{w) , u^g{w) are polynomial functionals and the function c is 
continuous. 

For E° the following bound is valid after adding to the argument of the exponent ^{w'g{P) — 
Wg{(3)) and applying the Schwartz inequality ( for the measure y/gdqP^gg^{dw)) 

e^" < glX, 

Jo = / exp{-2/3[(l - 3e)vg{w) + u.g{w)] + ^w'g\P)}Pfg{dw)dq, 

where wf{P) = J w'^{gt)dt. 


From the FK formula it follows that Jq is the trace of the kernel of the exponent of perturbed 
generator of the Wiener process. So, the Golden- Thompson inequality Tri^e"^^^) < Tr{e^e^) 
yields 

9I0 < V^(27r/?)-^ J exp{-2/?[(l - 3s)vg{q) + u,g{q) + U}dq = . (3.9) 

Here we took into account that 

exp{t9^}(g;g') = (477^)"^ exp{-^-^}. 

Iq is finite since h > Ad. 
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For /° after the rescaling q — (m ^g)^^q we have (g^ is the operator of multiphcation by g^) 



X exp{3"^(20eo ^ym-i^g^ + 8(m"^^)^ |g|)}c/g. 
Prom (3.3) it follows that 

exp{-^/3(--^a2 + ^g')}((m-^^)^g, (m-^^jig') = 
m Im Zg 

= e-^-^VI/^ ^M-^^i-d' + e - i)}(g, g') = 

= e-2-^^^exp{-^ + ^}exp{-yj/?(-ia2 + = 

= V^(l - exp{-^ + ^ - (1 - e-^^^)-\(^ - e'^^qf - 2-\f^P}. 

22 \ m 

Here we used the relation 

l{-d^ + f-l)^e-"^[~d^ + qd]e^ 

and the well-known formula for the density of the transition probability for the Ornstein- 
Uhlunbeck process. Hence 



X 



jO = ^^e-^\\ - e-2v^^)-5 / exp{3-^(20eo-^y^g' + 8(^)i|g|)} 
X cxp{-(l + e~V^^)"^(l - e''^)q^}dq. 

As a result 

7° = (1 - e-^v^^)-U(^)-i exp{f (3.10) 

y V TTX 

< ^e-v^(''-^'=(^)")(l - e-^^^)-h{g)-"n^, (3.11) 

where A;(5') is given by (1.9). 

Applying the Golden- Thompson inequality we obtain 

lu,, < V^Io- (3-12) 
Repeating (2.10-11), using the equality 

^ J Pf^idw) = V(27r/3)-im, 

we derive, also, the analogue of (2.12) for the quantum case 

{I^IX' ^ P^(irn-^r-'e^^^^^\ (3.13) 

As a result 

Combining all these bounds we see that e^{g) is bounded and (1.10) holds with 
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£;.(^)=lnJo- + e.(^). (3.14) 
Lemma 1.2 is proved. Theorem 1.1 is proved with an aid of Lemma 1.1 and (1.7). 



REMARKS. 

l.If one cancells the boundary term g then (1.1) is reduced to 

x&dA 

iSA \x—y\=l,x,y&A 

where dX is the boundary of X. If t/' > then the systems considered in [BF] can be recovered 
Surprisingly the proposed technique does not respect this boundary term since it creates an 
obstruction for obtaining the bound from above (2,6) for the rescaled and translated potential 
energy which guarantees uniform boundedness in g of lQu,g- 

2. If one cancells the boundary term 'U'iQxjQx) then (1.1) is equal to 

x€dX 

U{qx) = ^(^^' %) + U'{qx) u{q^, Qy) = {4d)~^{u{q^) + u{qy)) - gq^Qy 

x,yGX,\x—y\=l 

If U' is expressed in the same form as the first term in the rhs of the last equality then 
systems which are dealt with in [FL] can be recovered. Theorem (1.1) can be proved for such 
the potential energy taking into account in a special way a contribution of the boundary term 
to the sup erst ability, regularity conditions and (2.6) for the rescaled and translated potential 
energy (see [S]). 

3. Theorem 1.1 proves an existence of a phase transition for the case U' is expressed through 
a pair (special) potential, since it is known that in this case in the high-temperature phase there 
is an exponential decrease of correlations [K]. 

4. The magnitude of n-n interaction plays an exceptional role in the proposed approach 
since vanishing of it automatically implies vanishing of the spin two-point function for n-n sites. 
This means that E in (1.2) has to depend on the magnitude of n-n interaction, tending to zero 
together with it. So, one ought, always, to rescale by the magnitude (in an appropriate power) 
all the variables, when starting to derive the Peierls type contour bound using (1.7) with cq 
depending on it. 

5. Essentially ferromagnetic interaction may be characterized by the property that the fer- 
romagnetic configuration, consisting of the coordinate cq (minimum of a one-particle potential) 
at each lattice site, is more favorable than the associated antiferromagnetic (staggered) config- 
uration, consisting of the coordinate Cq at the even sublattice and — Cq at the odd sublattice 
for sufficiently large g, i.e. the potential energy on the former configuration is less than on the 
latter. This property follows from the superstability condition for the rescaled Ug in the formu- 
lation of Theorem 1.1 and the fact that the growth in g of g~^e^^ , s < n, is more slow than Cq. 
In other words, the ferromagnetic n-n part of the potential energy suppresses antiferromagnetic 
ground states for sufficiently large g. 

6. If one puts 

U\qK)= Cx-y{qx-qy)\ |C,_,| < ||C°||i<oo, 

x,y(iK 

where HC^Hi does not depend On g then the conditions of Theorem 1.1 are satisfied. 

7. The proof of Theorem 1.1 for classical systems and more general polynomial potentials 
u can be found in [S] 
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